Dwarf Spheroidal Galaxies in the Virgo Cluster 
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^ ; ABSTRACT 

! We present a study of the smallest and faintest galaxies found in a 

k> ■ ^^^y deep photographic R band survey of regions of the Virgo Cluster, to- 

u ■ tailing over 3 square degrees, made with the UK Schmidt Telescope. The 

■ " " ' objects we detect have the same physical sizes and surface brightnesses 

as Local Croup dwarf spheroidal galaxies. The luminosity function of 
these extremely low luminosity galaxies (down to Mr ~ —11 or about 
5 X 10~^L^) is very steep, with a power law slope a ~ —2, as would 
be expected in many theories of galaxy formation via hierarchical clus- 
tering, supporting previous observational evidence at somewhat higher 
luminosities in other clusters. 

Subject headings: galaxies: clusters: individual (Virgo) — galaxies: el- 
liptical and lenticular, cD — galaxies: evolution — galaxies: luminosity 
function, mass function — galaxies: photometry 
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1. Introduction 

Recent observations of rich clusters have 
indicated that the galaxy luminosity function 
(LF) may turn up at the faint end (e.g., Driver 
et al. 1994; Mohr ct al. 1996; Smith ct al. 
1997; Wilson et al. 1997; Trentham 1997). 
From a Schechter (1976) function slope of 
a ~ — 1, there appears to be a steepening to 
a. ~ —1.5 to -1.8 below about 0.04L^ (roughly 
Mr ~ -18 or Mb ^ -16.5 for Hq = 75 
kms^^Mpc^^). However the LF is still largely 
unknown faintwards of Mr ~ —15.5 {Mb — 
—14), so the overall contribution of dwarf 
galaxies to the cluster population remains un- 
certain. The Local Group provides the only 
well studied sample of such faint galaxies and 
it appears that below about Mr = —15.5 
(0.004L*) the galaxies are virtually all dwarf 
spheroidals (van den Bergh 1992a). These ob- 
jects span the magnitude range from about 
Mr ~ —14.5 (Fornax) to ~ —8.5 (Carina), 
and indicate a rather fiat luminosity distribu- 
tion (van den Bergh 1992b). 

In this paper we present a deep photo- 
graphic survey of significant areas within the 
Virgo Cluster which span a range of (gi- 
ant) galaxy densities. The survey, based 
on digitally stacking UK Schmidt Telescope 
films (cf. Bland-Hawthorn, Shopbell & Ma- 
lin 1993), extends the earlier, seminal, sur- 
vey of Binggeli, Sandagc & Tammann (1985), 
reaching roughly 3 magnitudes beyond their 
completeness limits. We are able to detect 
large numbers of faint galaxies, presumably 
dwarf spheroidals, at magnitudes down to 
Mr ~ — 11 for an assumed Virgo distance of 
18 Mpc (e.g., Jacoby et al. 1992; for consis- 
tency we also adopt their value of Hq = 75 
km s"^ Mpc~^ where required). 



2. The Data and Image Detection 

The data used here are part of a larger pho- 
tographic survey of the Virgo Cluster (Schw- 
artzenberg, Phillipps & Parker 1995a) using 
the extremely fine grained, highly efficient 
Tech Pan films on the 1.2m UK Schmidt Tele- 
scope (Phillipps & Parker 1993). Six indi- 
vidual long (1 to 1.5 hour) exposures of the 
same area, the South East quadrant of the 
Virgo Cluster, were scanned with the Super- 
COSMOS automatic measuring machine at 
the Royal Observatory Edinburgh (Miller et 
al. 1992). For convenience, nine separate 
scan regions 6840 pixels square were created 
from each film. The pixel scale is 10 mi- 
crons or 0.67" giving a total area for each scan 
~ 1.3° X 1.3°. Two of these scan regions are 
considered in the present paper, one near the 
cluster centre, one further out. Note that, 
while large in itself, the ~ 3.2 square degrees 
covered here is only about 1/40 of the full 
cluster survey area. 

The scans from the six separate films were 
sky subtracted by using a 256 x 256 pixel spa- 
tial median filtered version of the data them- 
selves, then matched in intensity by compar- 
ing images of a number of calibrating galax- 
ies and median stacked (see Schwartzenberg, 
Phillipps & Parker 1996 for details). Me- 
dian stacking has equivalent noise reduction 
to simple co-addition and is highly effective in 
removing artefacts (e.g., due to satellite trails, 
dust particles adhering to the emulsion and so 
forth) which affect only one film in the stack. 
Absolute calibration was via comparison of 
the images of some brighter galaxies with 
published CCD photometry (from Gallagher 
& Hunter 1989), as described by Phillipps & 
Parker (1993). The final stacked data have 
an equivalent exposure time of about 7 hours, 
and the high efficiency of the films, approach- 
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ing 10% (Parker et al. 1997), results in a 

pixel-to-pixel sky noise ag^y equivalent to 26.2 
R magnitudes per square arc second (hence- 
forth R^). 

Galaxy (and star) images were automat- 
ically recovered from the stacked data via 
a connected pixel algorithm (PISA, Draper 
1993), using a detection threshold 2a sky above 
the sky background, or 25.45i?/x, and a min- 
imum area of 25 pixels (11 square arc sec- 
onds). Each image thus has a minimum S/N 
of 10 and has a magnitude R <22. (In prin- 
ciple, 3(J detection of images of size around 
2" would reach a magnitude limit of about 
R = 24). Around 28,000 images are detected 
in each field and the large minimum area en- 
sures that few are spurious. This was con- 
firmed by comparison of a small area of the 
photographic data with a CCD image taken 
on the Anglo- Australian Telescope; 160 of 162 
images visible on the film were matched on 
the CCD frame, including all those in the 're- 
fined' samples used below. 

3. Results 

The two chosen for this study are 

centred close to M87 and 3.°1 to its south 
south east, enabling both the cluster core and 
a more typical cluster region to be investi- 
gated. The area immediately around M87 
was not searched due to its effectively higher 
background light level, so the inner and outer 
surveyed fields cover 1.58 and 1.61 square de- 
grees, respectively. Once the raw catalogues 
were produced, as above, we 'refined' them by 
requiring that our images met certain criteria 
aimed at isolating low surface brightness clus- 
ter dwarfs. In particular we kept only those 
images whose isophotal sizes and isophotal 
magnitudes were consistent with them having 
exponential profiles (characteristic of virtu- 



ally all dwarfs; see Binggeli & Cameron 1991) 
of scale size a > 2" and central surface bright- 
ness Hq > 22Rfi (cf. Figure 1 of Schwartzen- 
berg et al. 1995b). This reduced the num- 
ber of potential Virgo Cluster medium to low 
surface brightness dwarf galaxy candidates to 
approximately 17,000 (from 56,000 images of 
all types). Note that we do not separately 
remove stars as these should disappear along 
with the higher surface brightness galaxies. 
The data set includes galaxies with //q down 
to about 25.2 Rfi, but is complete (in the 
sense that even a 2" scale size gives images 
exceeding our area limit) only to 24.5i?yLi. Fur- 
ther details of other image parameters are 
given by Schwartzenberg (1996), but here we 
concentrate solely on the magnitudes, though 
note that we use 'total' magnitudes calculated 
from the measured a and fio- 

In principle it is possible for an LSBG sam- 
ple to contain cosmologically dimmed nor- 
mal surface brightness giants at large red- 
shifts or large non-cluster LSBGs in the back- 
ground. The former would generally appear 
much smaller than our detection limit (with 
a < 1", cf. Windhorst et al. 1994), while the 
latter are relatively rare (see Schwartzenberg 
et al. 1995b) and can be subtracted statisti- 
cally (Turner et al. 1993). Nevertheless, in 
order to reduce such contamination problems 
to a minimum, we have again refined our sam- 
ple to include only the ~ 4000 objects with 
a > 3". These images will also be less af- 
fected by seeing; even if the scale lengths are 
slightly increased by the blurring (and the rel- 
atively moderate resolution), the central sur- 
face brightness will be decreased to compen- 
sate, leading to little error in the derived total 
magnitudes. In effect we will have merely the 
LF of galaxies limited at a marginally smaller 
physical size than would have been the case 
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in the absence of seeing. 

Since we have pre-selected our dwarf LSBG 
sample in terms of scale size and central sur- 
face brightness, we cannot simply subtract 
standard number counts for the entire pop- 
ulation of field galaxies (e.g., Metcalfe et al. 

1995) from the magnitude distribution we ob- 
tain, in order to arrive at the cluster LF. We 
have therefore made a subtraction based on 
the corresponding distribution of field LSBGs 
parameters found by Schwartzenberg et al. 
(1995b). This correction turns out to be quite 
small compared to our total LSBG numbers 
(a few percent), so is not critical to our fi- 
nal LF, for the simple reason that most back- 
ground LSBGs appear much smaller than our 
cluster LSBG candidates (cf. Karachentsev 
et al. 1995). 

Figures la and lb illustrate the LFs for 
the outer (895 galaxies) and inner (675 galax- 
ies) cluster regions, respectively. We show 
here only those galaxies with a > 3" and 
22.0 < /iq < 24.5it!/i, the region of parame- 
ter space for which we have a complete and 
minimally contaminated sample. Note that 
this is not a magnitude limited sample. For 
instance, we already begin to lose any galax- 
ies with smaller scale sizes at Mr ~ —13, 
whereas our faintest objects have Mr ^ —11. 
Of course, even the loss of small objects at 
faint Mr may not be the whole story as far 
as the LF goes, since there may exist higher 
surface brightness dwarfs than we are allow- 
ing for (perhaps preferentially at bright Mr), 
and we will be missing any even lower surface 
brightness objects at all Mr. Indeed, recall 
that we have many candidates for smaller or 
lower surface brightness galaxies in our orig- 
inal overall sample (see also Schwartzenberg 

1996) . 

The LFs plotted for the two regions in Fig- 



ure 1 are very similar so an overall LF for the 
samples can be used. It is clear that the LF 
is again steep, as in the papers discussed in 
Section 1, confirming earlier suggestions for 
Virgo itself by Impey, Bothun & Malin (1988) 
and Tyson & Scalo (1988). A least squares 
fit to the combined data gives a power law 
slope for the range 15.5 < R < 20.0 (roughly 
-16 to -11.5 in Mr) oi a = -2.26 ± 0.13. 
(Fitting to the individual LFs over the same 
magnitude range gives « = — 2.26±0.14 and 
-2.18 ±0.12). Note that if we restrict at- 
tention to the very faint galaxies, i? > 18 
{Mr > —13.5), the steepening is even more 
dramatic, a ~ —2.5. (The turn up appears 
clearer in the 'outer' field. Figure la). The 
amplitudes of the LF for the two separate ar- 
also similar (in galaxies per magnitude 
bin per square degree), with the core sample 
actually having the lower projected density 
by a factor ~ 0.8. LSBGs may be adversely 
affected by the presence of the giant galaxy 
M87 in the cluster centre region; Thompson 
& Gregory (1993) have previously found a 
similar effect in the core of the Coma Clus- 
ter. Note, though, that with the very steep 
LF slope, a relatively small zero point offset 
in the calibration between fields, for instance, 
can have a significant effect on the numbers. 
For instance an error of Am — 0."*1 gener- 
ates a difference of a factor l0O-4(°+i)^"» ~ 
1.12. Such errors would make little differ- 
ence to the shape of the derived LF, since 
the background contamination is so small. 
At the bright end, the number of detected 
LSBGs is in good agreement with that ex- 
pected from the Binggeli, Sandagc & Tam- 
mann (1985, 1988) Virgo Chistcr LF, given 
the small numbers of these objects in our sam- 
ples (see Figure lb). It is clear that our LF de- 
parts from theirs at the expected point where 
incompleteness and lack of very low surface 
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brightness objects starts to affect tfieir sam- 
ple, beyond ^ 17. (We assume liere typical 
early type galaxy colours, B — R = 1.5, for the 
dwarf spheroidals; if they actually have bluer 
colours, as often seen in low surface bright- 
ness galaxies, this would slightly improve the 
match) . 

4. Conclusions 

By co-adding very deep UKST photographic 
films we have been able to probe the dwarf 
population of the Virgo Cluster down to Mr ~ 
-11 (~ 5 X 10~^L^). The central surface 
brightness limit for our sample is 25i?/i, cor- 
responding roughly to 26.5Bfi for early type 
galaxy colours. In both luminosity and sur- 
face brightness this is thus one of the deep- 
est surveys yet performed. In particular, 
our limits allow us to survey well into the 
regime of the dwarf spheroidal galaxies (Irwin 
& Hatzidimitriou 1995); the luminosity limit 
is 25 times fainter than the Fornax dwarf, 
for instance. We have therefore been able to 
gather by far the largest sample of dwarf ellip- 
ticals/dwarf spheroidals currently known. Of 
course, in the absence of redshifts, these 'de- 
tections' are on a statistical basis only. How- 
ever, given the paucity of LSBs of moderate 
to large angular size in the general field, we 
are confident that the large majority of our 
candidates are genuine cluster dwarfs. 

The luminosity function of the dwarfs is 
very steep, with a ~ —2, confirming val- 
ues found over much more limited magni- 
tude ranges in other clusters (e.g.. Smith et 
al. 1997; Trentham 1997). Bernstein et al. 
(1995) reached similarly faint levels to those 
discussed here with very deep CCD imaging of 
a very small area at the core of the Coma clus- 
ter (to Mr — —11.4). They found a less steep 
LF than most other deep surveys, a ~ —1.3, 



but the centre of Coma may be a rather spe- 
cial environment. Though their surveys are 
less deep, Biviano et al. (1995) and Thomp- 
son & Gregory (1993) find steeper slopes for 
parts of Coma further from the centre. 

A steep slope, a ~ —2, is as expected gen- 
erally in any hierarchical structure formation 
model (eg. White & Frenk 1991; Blanchard, 
Valls-Gabaud & Mamon 1992; Evrard, Sum- 
mers & Davis 1994; Frenk et al. 1996; Kauff- 
mann Nusser & Steinmetz 1997). Note that 
since we are observing at long wavelengths 
{R band), and in any case we expect most 
of our objects to be dwarf ellipticals with lit- 
tle or no recent star formation, our LF shape 
should closely match that of the more funda- 
mental (baryonic) mass function, allowing a 
simpler comparison with theoretical models. 
There is a suggestion that the dwarf LSBG 
to giant galaxy ratio is smaller in the cluster 
core than further out. Analysis of the whole 
cluster survey area should allow us to quan- 
tify this in more detail (see also Phillipps et 
al. 1997). 

We might, finally, note that much smaller 
and fainter galaxies can be detected in our 
data than are present in our photometric sam- 
ples. Indeed, we can reach down to about 
Mr ~ —8.5 at the distance of Virgo, the 
same as the Carina dwarf, the lowest luminos- 
ity system currently known. Unfortunately, 
though, these images are indistinguishable 
from those of the (very numerous) general 
background galaxies. However, we should still 
be able to estimate their numbers through a 
comparison with an identically observed non- 
cluster field. This work will be reported in a 
subsequent paper. At some point we should 
certainly expect to see a turn over in the LF, 
since a = —2 is the critical value at which 
the integrated galaxy light formally diverges. 
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If the currently found slopes in the range -2 to 

-2.5 were to continue down to, say, Mr = —8, 
then the dwarf galaxies fainter than = 16 
would contain approximately 0.1 to 1.0 times 
as much hght as the brighter galaxies. For a 
constant M/L (or perhaps more reasonably, 
a fixed baryonic M/L) this would obviously 
increase the total mass in cluster galaxies by 
a factor between 1.1 and 2. 

We would like to thank the UK Schmidt 
Telescope for the provision of the usual ex- 
cellent photographic material and the Super- 
COSMOS group at ROE for scanning them. 
SP and JBJ thank the Royal Society and the 
UK PPARC, respectively, for financial sup- 
port. 



REFERENCES 

Bernstein G.M., Nichol R.C., Tyson J.A., Ul- 
mer M.P., Wittman D., 1995, AJ, 110, 
1507 

Binggeh B., Cameron L.M., 1991, A&A, 252, 
27 

Binggeli B., Sandage A., Tammann G.A., 
1985, AJ, 90, 1681 

Binggeli B., Sandage A., Tammann G.A., 
1988, ARA&A, 26, 509 

Biviano A., Durret F., Gerbal D., Le Fevre 
O., Lobo C, Mazure A., Slezak E., 1995, 
A&A, 297, 610 

Blanchard A., Valls-Gabaud D., Mamon 
G.A., 1992, A&A, 264, 365 

Bland-Hawthorn J., Shopbell P.L., Malin 
D.F., 1993, AJ, 106, 2154 



Driver S.P., Phillipps S., Davies J. I., Morgan 
I., Disney M.J., 1994, MNRAS, 268, 393 

Draper P, 1993, PISA Manual, STARLINK 
User Note No. 109 

Evrard A.E., Summers F.J., Davis M., 1994, 
ApJ, 422, 11 

Frenk C.S., Evrard A.E., White S.D.M., Sum- 
mers F.J., 1996, ApJ, 472, 460 

Gallagher J.S., Hunter D.A., 1989, AJ, 98, 
806 

Impey C, Bothun C, Malin D., 1988, ApJ, 
330, 634 

Irwin M.J., Hatzidimitriou D., 1995, MN- 
RAS, 277, 1354 

Jacoby G.H., et al., 1992, PASP, 104, 599 

Karachentsev I.D., Karachentseva V.E., 
Richter CM., Vennik J.A., 1995, A&A, 
296, 643 

Kauffmann C, Nusser A., Steinmetz M., 
1997, MNRAS, 286, 795 

Metcalfe N., Shanks T., Fong R., Roche N., 

1995, MNRAS, 273, 257 

Miller L.A., Cormack W., Paterson M., Beard 
S., Lawrence L., 1992, in MacGillivray 
H.T., Thomson E.B., eds. Digital Optical 
Sky Surveys, Kluwer, Dordrecht, p. 133 

Mohr J.J., Geller M.J., Fabricant D.G., Weg- 
ner G., Thorstensen J., Richstone D.O., 

1996, ApJ, 470, 724 

Parker Q.A., Phillipps S., Malin D.F., Can- 
non R.C., Russell K.S., 1997, MNRAS, to 
be submitted 



6 



Phillipps S., Driver S.P., Couch W.J., Smith 
R.M., 1997, ApJ, to be submitted 

PhiUipps S., Parker Q.A., 1993, MNRAS, 265, 
385 

Schechter P., 1976, ApJ, 203, 297 

Schwartzcnbcrg J.M., 1996, Ph.D. Thesis, 
University of Bristol 

Schwartzcnbcrg J.M., PhiUipps S., Parker 
Q.A., 1995a, A&A, 293, 332 

Schwartzcnbcrg J.M., Philhpps S., Parker 
Q.A., 1996, A&AS, 117, 179 

Schwartzenberg J.M., PhiUipps S., Smith 
R.M., Couch W.J., Boyle B.J., 1995b, MN- 
RAS, 275, 121 

Smith R.M., Driver S.P., Phillipps S., 1997, 
MNRAS, 287, 415 

Thompson L.A., Gregory S.A, 1993, AJ, 106, 
2197 

Trentham N., 1997, MNRAS, 290, 334 

Turner J. A., Phillipps S., Davics J. I., Disney 
M.J., 1993, MNRAS, 261, 39 

Tyson N.D., Scalo J.M., 1988, ApJ, 329, 618 

van den Bergh S., 1992a, MNRAS, 255, 29p 

van den Bergh S., 1992b, A&A, 264, 75 

White S.D.M., Frenk C.S., 1991, ApJ, 379, 52 

Wilson G., SmaU I., Ellis R.S., Couch W.J., 
1997, MNRAS, 284, 915 

Windhorst R. et al, 1994, AJ, 107, 930 



This 2-column preprint was prepared with the AAS 
IMgX macros v4.0. 

7 



16 



14 



12 



1000 



100 



10 



10 



~| \ \ 1 \ 1 \ 1 1 \ 1 1 1 r 

- a = -2.26±0.14 



_l I L_ 



1000 r 



CO 



OB 
CO 

:.0 



20 



22 



100 



10 r 



16 
1 ' 



14 



12 



10 



~i \ r 



a = -2.18±0.12 



J I L_ 



20 



22 



Fig. 1. — The absolute magnitude distribution of Virgo LSBGs (as defined in the text), for the 
(a) outer and (b) inner area samples. Error bars shown are based on Poissonian statistics. The 
luminosity function of Sandage et al. (1985) is also shown (open triangles) for the Virgo Cluster 
Catalog galaxies which overlap with the inner field. (We assume for this comparison B — R = 1.5.) 
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